clear all; close all;

% delt = 1000/179857.92;
delt = 1e-4;
N = 5000;
t = (0:N-1)*delt;
W = -0.02 + 0i;
s = exp(W*t);
WI = imag(W);

L = 1000;
z = zeros(1,L);
x = zeros(1,L);
for ii = 1:L
    ww = real(W) - 50i + ii*.1*i;
    x(ii) = imag(ww);
    A = s - exp(ww*t);
    z(ii) = A*A';
end
plot(x,z);